3D two-color QCD at finite temperature and baryon density 
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We study the phase diagram for two-color QCD in three-dimensional spacetime, as a function of 
temperature and baryon chemical potential, using the low-energy effective Lagrangian approach. We 
show one-loop renormalizability at zero temperature, and then use the one-loop effective Lagrangian 
at finite temperature and chemical potential to show that at low temperature there is a critical line 
separating the normal and diquark phase, with this critical line ending at a tricritical point. This 
phase structure is qualitatively similar to that found recently by Splittorff et al. for two-color QCD in 
four-dimensional spacetime, although the details are quite different, due to the different symmetries 
and the different loop and infrared properties of three-dimensional spacetime. 

PACS numbers: ll.10.Kk, 11.30.Hv, 11.30.Qc, 12.39.Fe 

Keywords: chiral Lagrangian, chemical potential, finite temperature, baryon superconductivity 

I. INTRODUCTION 

A systematic approach to the investigation of the symmetry breaking structure of fermion-gauge theories is to 
use low-energy effective Lagrangians 0, to study the Goldstone modes corresponding to the spontaneously broken 
v—i • global symmetries |3] . Ideally, one would like to understand the phase diagram of fermion-gauge theories in the (T, fx) 
^ \ plane, where T is the temperature, and \i is the (baryon number) chemical potential. However, it is not known how 
to analyze general fermion-gauge systems at finite baryon density, since the baryon number chemical potential makes 
the Euclidean Dirac operator non-Hermitian and the Boltzmann weight complex. This problem can be avoided, as 
advocated in 0, H, 0,0 @ , by considering "two-color QCD", for which the fundamental representation of SU(2) is 
pseudoreal. From studies of 4D theories, it has long been appreciated that the SU(2) theory exhibits exotic types 
£f*^ , of spontaneous breakdown of global symmetry 0. Quarks and charge-conjugated antiquarks are combined into an 
extended flavor multiplet, which is expected to break into its extended vector subgroup. This has the consequence 
that in lattice simulations the Boltzmann weight is real even at finite (baryon number) chemical potential ll. For an 
even number of flavors, the Boltzmann weight is also positive definite. Thus, analytic predictions can be quantitatively 
compared with Monte Carlo simulations in lattice gauge theory [Tol [TT| , provided the lattice regularization respects 
Oh, the relevant flavor symmetry group. Such studies have recently been carried out |12| in 4D QCD with quarks in 
pseudoreal (and real) representations at finite fx. 

In ll| we studied the T = case of 3D two-color QCD, and found that the phase structure predicted by the 
tree-level effective potential was qualitatively the same as the 4D case studied in y| la, 13 , even though the details of 
the symmetry breaking in 3D are completely different from in 4D. This universality between 3D and 4D corresponds 
to the existence of two different breakings of the original global symmetry by a mass and a chemical potential term. In 
this paper we extend the 3D analysis of [13| to include nonzero temperature effects via the one-loop effective potential, 
as was done recently for the 4D system in [l^. The loop effects introduce another difference between 3D and 4D, 
as infrared physics is more significant in 3D. As a result, we find some similarities with the 4D case, but also some 
important differences. 

In 3D, with an even number Np of flavors of massless complex fermions, denoted by the Np/2 pairs ipt, Xf, one can 
predict spontaneous flavor symmetry breaking patterns along similar lines to 4D QCD . In 3D the generic situation is 
that the continuous part of the global symmetry group is broken according to [la , Il7| 



43 
Oh 



43 



- 1—1 

X 

5-H 



V(N F ) -> U(N F /2) x U(AT F /2) (1) 



by the quark-antiquark condensate 



N F /2 

E ((Mf) - (XfXf)) , (2) 
/=i 

Evidence for such a symmetry breaking pattern has been observed in 3D lattice simulations 01 with gauge group 
SU(3). This pattern of flavor symmetry breaking can also be predicted for 3D QCD at large Nc using the Coleman- 
Witten argument [l9|. The 3D symmetry breaking pattern in (JIJ is for quarks in a complex representation of the 
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gauge group, and is expected to apply to a U(l) theory and to SU(iVc) theories with number of colors Nc > 3. 
For SU(2), with fundamental quarks, the symmetry breaking pattern is expected to be different again, due to the 
pseudoreality of the fundamental representation. The pseudoreality of the fundamental representation of SU(2) means 
that the naive \](Np) flavor symmetry is extended to USp(2A^), and the continuous part of this global symmetry 
group is predicted to break down in 3D as poj 

USp(2V F ) _> USp(A^) x USp(V F ). (3) 

This is different from the generic 3D symmetry breaking pattern in JQ), and also is different from 4D theories where 
the standard flavor symmetry breaking patterns are 

V(N F ) L x \J(N f )r SU{N F )v , for SU(V C > 3) 

U(2V F ) -> USp(2V F ) , for SU(2) . (4) 

In 4D the U(iV F ) flavor symmetries are first broken to SU(iV F ) symmetries by the axial anomaly, and then broken by 
the chiral condensate, with the net breakings as shown in (J2J. Physically, the differences between the 3D and 4D cases 
reflect the differences between the anomalous discrete symmetries of parity and chirality in 3D and 4D, respectively. 

In section II we review the global symmetries and their cosets that are relevant for the low energy effective Lagrangian 
description of the 3D two-color QCD system, with quarks in the fundamental representation. In section III we compute 
the one-loop effective action at finite temperature and finite baryon chemical potential. This result is used in section IV 
for a Landau- Ginzburg description of the phase diagram, using the rotation angle a which characterizes the tree-level 
vacuum state as an order parameter. Section V contains our conclusions and some discussion. 



II. SYMMETRIES AND EFFECTIVE LAGRANGIAN 

A. Enlarged flavor symmetry 

The fermionic part of the microscopic Lagrangian for two-color QCD with Np = 2n flavors of quarks in Euclidean 
3D space is 



C f = ifjftip + xflx + mitpip - xx) ~ MW'V + X f x)- 



(5) 



Here ip = ip^, if)' = ip* l f, X = X/? X^ — X*/ are independent two-component spinor fields, with the color index i = 1, 2 
and flavor index / = 1, ...,n being suppressed. Pauli matrices denoted as a v , with v = 1,2,3, are employed to 
represent the Euclidean Dirac matrix algebra, and those denoted as r Q , with a — 1,2,3, are employed to represent 
the gauge group algebra. We choose X3 to be the Euclidean time direction and define ip = ip^&s, x — X CT 3- The Dirac 
operator is If) — cr v D U: and the covariant derivative is D y = d v + iA V: and the gauge field A v = A^r a is Hermitian 
and su(2) valued. 

Due to the pseudoreality of the SU(2) Dirac operator, the kinetic part of the Lagrangian © is invariant under a 
symmetry group larger than the apparent U(2n) symmetry. One can put i/j, x, tp, and x into a single flavor 4n-plet, 



f = 



1> 
X _ 

CT2T2lp T 

a 2 T 2 x T 



(6) 



so that 



(7) 



where the An x 4n matrices /, M and C are 
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-t 2 „ ' 


_ t 2n 







-In 




In 






-In 





c = 





t 2n 


t 2n 





(8) 



3 



The kinetic term in the microscopic Lagrangian (JJJ has extended flavor symmetry group being the unitary symplectic 
group USp(4n), 



S T IS = I, S^S = 1 4 



(9) 



The above extension of the flavor symmetry group is analogous to what happens in two-color QCD in 4D, where 
the conjugated right-handed spinor <J2Ti' l l } \ transforms as the left-handed spinor ipi does under gauge and Lorentz 
transformations, so that the chiral S\J(Nf)l x SU(Nf)r symmetry is extended to SU(2Np). 

In the 3D case, the flavor group USp(4n) is broken down to USp(2n) x USp(2n) explicitly by the mass term in J7J, 
or spontaneously by the quark-antiquark condensate J2J, if formed. In the latter case, the Goldstone manifold is thus 
given by the coset space USp(4n)/(USp(2n) x USp(2n)), which has An 2 independent degrees of freedom, and can be 
parametrized as 



where 



S(x) — exp 



ill(x) 
2F 



Ii{x) = TT a (x)X a . 



(10) 
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The fields ir a are the Goldstone modes, and the 4n 2 generators X a span the subspace usp(4n) — (usp(2n) © usp(2n)). 
The choice of E c = M> leads to the block representation of the generators X a as 



(12) 



where P and Q are n x n complex matrices, each having 2n 2 degrees of freedom. The chemical potential term in (|7|) 
explicitly breaks the extended flavor group USp(4n) down to its unextended U(2n) subgroup, while in the presence of 
both mass and chemical potential terms, the surviving global symmetry becomes U(n) x U(n), which is the intersection 
of USp(2n) x USp(2n) and U(2n). 



B. Low-energy effective Lagrangian 



Coupling the fermionic Lagrangian Q) to the SU(2) gauge field gives the the microscopic Lagrangian 

1 



(13) 



In the low energy Aqcd 3 ) regime, where fundamental particles are confined and Goldstone bosons dominate, we 
can define a low energy effective theory. The kinetic term of the effective Lagrangian describing the Goldstone modes 
E, parametrized as in I|1U|) and Ijllfl . should be invariant under the action of the global USp(4n) group 



-> sS(x)s T , s e USp(4n) 



(14) 



with s in the antisymmetric tensor representation. Invoking the standard symmetry principles of chiral perturbation 
theory, and the principle of local flavor symmetry, the effective Lagrangian is 



L = — tr V^SV^E 1 - F 2 M 2 tr (ME) 

where the flavor-covariant derivatives are 

V y E = d v Il - n(Y.Bl + BlT) 
V,E t = djtf + B v + B v Yt). 

with 



B = CI 



l 2 n 






~t 2 n 



(15) 



(16) 



(17) 
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C. Phases at tree level 

At zero temperature, the static part of the effective Lagrangian (JTBJ, i.e. the effective potential, determines the 
vacuum alignment of the field S. Introducing a dimensionless parameter £ = 2fx/M to represent the chemical potential, 
we find 



L st (S) = F 2 M 2 (- tr (MS) - L- tr (BY,Btf) - 



(18) 



The above two terms in S compete for the direction of the condensate which we denote by S. The condensate that 
gives the global minimum of the static tree level effective Lagrangian is 

S Q = cos a 4- J sin a, (19) 
cosa = min(l,£~ 2 ). (20) 



which is parametrized by a where 



1. Normal phase 

When £ < 1, (i.e., /x < 4f ), the vacuum orientation of the condensate does not depend on £ and is given by S = M< . 
Expanding £ around AP using the Goldstone field defined in (|1 If) and (|12|) we find 

L = L rt (£ a ) + - tr d v P^B v P + — tr P*P + - tr d v Q?d v Q - 2^ tr Q%Q + - 2/z 2 J trQ t Q + ---. (21) 

From this we can identify the excitations. The baryon charge b of the excitations are : b — for P and P>; 6 = 2 for 
Q; and b = — 2 for Qt. The pole masses for these excitations are 

M for P and P f 

M-2/i for Q 

M + 2/i for Q f . (22) 
Each mode belongs to a representation of dimension n 2 . 

£. Diquark condensation phase 

In the regime £ > 1, (i.e., /i > 4f ), the vacuum condensate has U(n) degeneracy corresponding to n 2 true Gold- 
stone modes. This change of massless modes indicates a second-order phase transition at £ = 1. When £ > 1 the 
configuration Ijl9(l begins to rotate away from ilP according to cos a = £~ 2 . This rotation can be written as 

t a = s a MU T a = s 2 a M\ Sa = e^ x \ (23) 

where X2 is the generator that rotates M into J . We parametrize the fluctuation around the vacuum £„ as 

S = SaSM^SFsl, (24) 

where S 1 are generated by the unrotated generators X a . 

In this phase we expand around the rotated value of the condensate £ a , as in H24J) . To second order 

L = L 8t (E ) + ~ tr dvPldvPs + \ tr d v P\d v P A + X - tr (^Q fl ) 2 + 1 tr [d v Qif 

M 2 r i 
+ iMr 1 tr(Q^ 3 Q/-Qla 3 Q«) + — [£ 2 trP s pt +r 2 tr p A p^ + ^2_ r 2 )trQ2 j + ..._ (25) 

Here Ps^ = (P ± P T )/2, and Qa,/ are the real and imaginary parts of the complex matrix Q. The Ps,a fields have 
conventional dispersion relations with masses 

2fi for P s 

M 2 

■Tj- for Pa (26) 
2/x 
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and lie in representations of dimension n(n + 1) and n(n— 1), respectively. 

On the other hand, the dispersion laws for the Q fields are unconventional due to the linear time derivative terms 
in and are determined by the secular equation 



det 



[E 2 - p 2 ] 2iEM£- 1 

-21EMC 1 [£ 2 -p 2 -M 2 (e 2 -r 2 )] 



(27) 



Diagonalizing, we find one true Goldstone field, which we denote Q, and its orthogonal complement denoted . 
These excitations have masses 

for Q 



M N 4 



Wl + 3^j for & (28) 

and each lies in a representation of dimension n 2 . Note that fields with this unconventional form of dispersion also 
appear in the context of kaon condensation 0,H2, and gauge symmetry breaking via Bose condensation |23|. 

III. ONE-LOOP FREE ENERGY AT FINITE DENSITY AND TEMPERATURE 

The bosonic low energy effective Lagrangian for 3D two-color QCD reads 

F 2 

L = —trW„Y,(x)W^{x) i -F 2 M 2 trMJ:(x)+L 4 . (29) 

where the phenomenological constant F 2 has dimensions of mass in 3D, and L4 is a collection of terms of order 
O (V, M) 4 . We take M > without loss of generality. 

A. Normal phase 

Within the normal phase (the boundary of which will be determined in the next section) , the Lagrangian reads 

L = -AnF 2 M 2 + ip(a;)t(_a 2 + M 2 )P{x) + ^Q(x^(-d 2 - 4/i«9 3 + M 2 - V)Q(x) + L A . (30) 

Thus the one-loop free energy is given by 

Q = -AnF 2 M 2 + n 2 tr log(-9 2 + M 2 ) + n 2 tr log(-<9 2 - Afid 3 + M 2 - 4/i 2 ) + L 4 [S = M f ] 

= -AnF 2 M 2 + n 2 ^^log(p 2 + M 2 ) + n 2 \og(p 2 + U m + M 2 - V) + Sh- (31) 

Here LI4 is a constant that is a combination of the contact term couplings. The sum-integral in the above is defined 
as 

r. OO p 

Td d pf( P3 ,p) = 2TrT £ d d - 1 pf(2nnT,p). (32) 

n=-oo^ 

The first sum-integral in 131|l is standard 



d=3 M 3 , T 
' 6 



' N=l ^ ' 

- + — fMLi 2 (e-^)+TLi 3 (e-T)). (33) 

7T 7T \ / 
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Note that the zero temperature part is finite in odd dimensional spacetime, and the finite temperature correction is 
expressed in terms of the polylogarithm function [24[ 



U n (x) = 



k=l 



(34) 



The second sum-integral in (|31|l is the conventional effective action for a bosonic field at finite chemical potential, 

d d p 



t 



(2tt)< 



■log((p 3 + 2inf 



1 



( 47r )d/2 



dt 



I E 



M 



e 4tT 2 exp 



\ T J 



(47r) d / 2 

=3 M 3 T 2 r 



V 27 



d/2 



N=l 



MN 
T 



cosh 



d=3 



6tt 2-7T 



M 



f\U 2 (e — ) + Li 2 (c —)) +T(Li 3 (c — ) + Li 3 (e — )J 



(35) 
(36) 



Once again, the zero T part is finite due to the odd dimensionality, and the finite T corrections are expressed in terms 
of polylog functions. 

Bringing together the contributions (|33|l from the P modes, and the contribution Ij36(l from the Q modes, the 
onedoop effective action in the normal phase, at finite T and fj,, is 



n 



-4nF 2 M 2 



M 3 

67T 



— (MLi 2 (e--)+TLi 3 (e-- 

7T V 



(M 3 


T 2 








M ( 


\ 6tt " 







M-2 M M + 2 M 

'^^) + Li 2 (e ^ 



, Af-2 M , M + 2 M , 

T Li 3 (e ^)+ Li 3 (e ^) 



(37) 



B. Diquark condensation phase 



We now compute the free energy, to ldoop order, by expanding about the rotated vacuum 

cxp Of X2 



(38) 



which minimizes the tree level effective potential in the diquark condensation phase. We parameterize the fluctuations 
around E Q as 

E(x) = s a S(x)M^S(x) T sl, 

'm(x 



S(x) 



exp 



2F 



, n(a:) 





P 




Q 






Q T 






Q* 










-P T 





(39) 



The quadratic Lagrangian reads 

L = -F 2 fi 2 ti[BE a BEl + Un) - F 2 M 2 tiMH a 



+ l -P s {x)\~d 2 + Ml + \m 2 3 )P s {x) + l -P A { x )^-d 2 + M 2 2 + \m 2 3 )P a {x) 



1 



-d 2 + Ml M 3 d 3 
M 3 d 3 -d 2 + Ml 



Q{x) +L 4 . 



(40) 



Here we have used the notations of expressing Mi, M 2 and M 3 in terms of the overall mass scale M, the 

baryon chemical potential /Lt, and the rotation angle a which characterizes the rotated vacuum (|38|) : 



Ml = M 2 cos a - 4[i 2 cos 2a, 
Ml = M 2 cos a - 4^ 2 cos 2 a, 
Ml = 16a/ 2 cos 2 a. 



(41) 
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There are n(n + 1) P$ modes, n(n — 1) Pa modes, and (n 2 + n 2 ) Q modes. Integrating out these modes, one obtains 
the one-loop free energy 



CI 



F V tr {BZaB^i + t in ) - F M tr MS a 

X -n(n + 1) tr log(-d 2 + M 2 + ^M 3 2 ) + \n{n - 1) tr log(-d 2 + M 2 + Jm 3 2 ) 



1 



-n tr log 



-d 2 + M 2 



M 3 d 3 



M 3 8 3 -d 2 + M 2 



-AnF 



2 (M 2 + M 2 Mi 



i(n — 1) srjk d 3 p 



n(n + 1) v-A d 3 p . , „ , 1 „ ~. 



1 1 M\ 



-Mi 



(42) 



Here 514(a) = LiX^ a ] is given by substituting T,(x) — S Q into the 0(V,M) terms of the effective Lagrangian. 

The Ps and Pa mode contributions in (|42|) are precisely as for the P modes in the normal phase in (|33|) . with M 
replaced by \J M 2 + M 2 / 4 for the Ps modes, and by y/ M% + Mf / 4 for the modes. On the other hand, the final 
term in l|42|l is nontrivial because of the quartic nature of the operator. One approach is to factor this quartic as 



(p 2 + M 2 ){p 2 + Mi) + P 2 3 M 2 = ((p a + iMy) 2 + p 2 + M 2 z 2 ){(p 3 - iMy) 2 + p 2 + M 2 z 2 )) 
where we have introduced the convenient dimensionless parameters y and z, as in [bH Il5j. 



(M 2 -M 2 



(43) 



1 / Ml + Mi M 2 



7 = ^ 

V 2M ' Z M V 2 -I 
One can now expand the final term in (|42|l in terms of the difference 

Ml - Ml = Ap? sin 2 a. 



(44) 



(45) 



Such an expansion proves useful in the Landau- Ginzburg analysis for the region where the order parameter a is small, 
as studied in the next Section. 

Then the contribution of these Q-type modes is 



Gn = — 



log {{p 2 + M 2 ){p 2 + Ml) +p 2 M 2 ) 

log (((p 3 + zMy) 2 + p 2 + M 2 z 2 )((p 3 - iM y) 2 + p 2 + M 2 z 2 )) 



4^(2tt)° 



Y- 

fc=i 



1 (Ml -Ml 



2 k 



) d ((p 3 + iMy) 2 + p 2 + M 2 z 2 ) k ({p 3 - iMy) 2 + p 2 + M 2 z 2 ) k ' 



(46) 



The logarithmic term in the second line of (|46|l factors in terms of conventional bosonic effective actions at finite 
chemical potential : 



G 



(0) 



log [((ps + iMy) 2 + p 2 + M 2 z 2 )((p 3 - lM y) 2 + p 2 + M 2 z 2 )} 



dt- 



-t(Mz 



(47r) d / 2 J Q t d / 2 + 



— e itT* cosh 



N=-c 



MyN 
T 



T(-d/2), Nrf /MzT\ d/2 ^ rf/ „ (MzN\ , /MwiV 



(47) 



(2=3 



(Mz 



3tt 



|3 y2 
7T 



Mz(Li 2 (e- M ^)/ T ) + Li 2 (e- M ^)/ T )) + T ( Li 3 ( e - M ^)/ T ) + Li 3 (e~ M ( z+ ^/ T )) 
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The finite-temperature parts of l|47[) correspond to the corrections for masses M(z ± y). 

The k > 1 terms in the expansion (|46[) are more difficult to evaluate in a simple closed form. A natural approach 
is to use Feynman parameters: 

vf d d P 1 

4^ (27r) d ((p 3 + iMy) 2 + p 2 + M 2 z 2 ) k ((p 3 - iMy) 2 + p 2 + M 2 z 2 ) k 



r(2fc) vfv d d p Z" 1 



(*(!-*)) 



fe-i 



T(fc) 2 4^ (27r) d 7 {p 2 + M 2 (z 2 - y 2 ) + (2s - l) 2 iMyp 3 ) 2k 

-| pi poo x 

^—j^ / ds^l-a))*- 1 / dti 2 *-i- i! re- tM3 (* > -» a )T V 

(ATT) 2 Jo JO r, — r 



-t[{2irnT) 2 + (2s-lf iMy(2irnT)} 



n— — oo 



r(fc) 2 (47r)- 

1 pi rOO 00 9 

— -1— - / ds(s(l- S )) k - X dU 2k-l- d /2 e -tM^-y>) e -tM^ { l-2,f y ^^ + ^(1-2. 



r(fe) 



2 jfe V™/ Jo 



' z 2 — y 2 u 



d=3 2 2 4K 7T 2 



E 



2 A:- 



1 



,fc-i Jf £ 



4-2fc 



(z 2 — y 2 u) k 4 



cosh 



AT=— OO 

f MN 



z 2 — y 2 u 



(z 2 — y 2 u) k J 



cosh 



V T 
{MN 



yVT 



(T = part) 



V 71 



yVl^u + (T = part) 



(48) 

Motivated by the Landau- Ginzburg analysis of the next Section, we consider just the k = 1 term, since higher terms 
are at least of order a 8 . In d — 3 the fc = 1 term in (|48|l reduces to a simple integral : 



8irMy 1 2 2 - y J( z - v )m/t x {e x ~ 1) 



(49) 



Here the first term corresponds to the T = contribution, while the integral represents the finite T correction. 

Collecting together the contributions of the Ps, Pa and Q modes, we obtain the following expression for the one-loop 
effective action in the diquark phase, at finite temperature and chemical potential 



n = -4nF 2 M 2 z 2 



n(n+ 1) l" (M 2 + Af 2 /4) 3 / 2 T 2 
2 6n 7r 

n{n-l) (Mf + Mf/4) 3 / 2 T 



Ml + ^ U 2 (e^ M ? +M "^ T ) + TLi 3 (e-^ 1 2 +^3 2 /4/T ) 



6tt 



^ ^M 2 + M Li 2 (e-V M l+A^/4/T ) + rLi3(e -VM|+M|/4/T^ 



n 2 r m 3 z 3 r 1 

(M 2 - M 



Mz ( U 2 (e-^ M ' T ) + U 2 (e-^ M ' T )\ + T (Li 3 (e-^ M / T ) + Li 3 (e-< Z+ ^ M / T ) 

+ o(m\ - m|) 4 + n 4 . 



2 l0g — + J iz - y)M/T i^l) 



327ryM 

Note that the higher order corrections are at least of order (M 2 — M|) 4 ~ /x 8 sin 8 a 



(50) 



IV. LANDAU-GINZBURG MODEL 



In the Landau- Ginzburg approach |25|. we regard the rotation angle a as an order parameter. For small a, the 
shape of the effective potential can be used to deduce information about the phase structure of the system [26| . The 
vanishing of the coefficient of the term quadratic in a determines a critical line in the T and /i plane. 



A. Landau-Ginzburg model at zero temperature 



It is important first to analyze the Landau-Ginzburg model at T = in order to understand the effect of one-loop 
renormalization. Since we expect the critical behavior to be near the tree-level T = critical point where [i = M/2, 
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we expand f2 in terms of the rotation angle a and also the deviation of chemical potential from M/2: 

Substituting the definitions of Mi, M2, M3, z and y into the T = parts of the one- loop effective potential in 15(J|) 
and expanding, one finds 



= est. + -±- " a» + 2 + W + O (a 4 log(a 2 - 40) . (52) 



^t-o 1 M o / 2n + lM \ , 

~ = = est. H ^a 2 + 2 H ^ > --'" 2 

—AnF 2 M 2 32tt F 2 \ 87r F 

Note that the higher order terms in l|50(l are not relevant as they contribute terms at least of order a 8 . Also, we 
neglect the quartic contact terms Q4. This is motivated by the 4D case [l^ where these terms are negligibly small. 
The dimensionless combination M/F 2 in (|52|) is the loop expansion parameter. The vanishing of the coefficient of a 2 
leads to the critical point (at one-loop) 



M = Ml-iH) (54) 



This can be rewritten as 

2^ V" 32tt F 2 J 

which shows that the critical point is shifted by a finite one-loop correction from the tree-level T — critical point 
fi = M/2. Actually, this shift is just a mass renormalization since the one-loop renormalized pion mass m n is 



8(47r) d / 2 F 2 J \ 32irF 2 J 

as computed in the Appendix. Thus, this result means that one-loop renormalization has the effect of shifting the 
T = critical point from fj, = M/2 to (i = m n /2, where m v is the one-loop renormalized value of M. 



B. Landau-Ginzburg model at finite temperature 

At finite temperature, the shape of the free energy as a function of the order parameter a changes as we move 
through the T and \x plane. Given the explicit expression ifBTTjl it is straightforward to plot (using, e.g., Mathematica) 
the free energy for chosen values of T and /1. This gives some guidance as to the parameter region in which to search 
for a critical line. We find that the appropriate critical region of parameters is 

< (i ~ m n g <C T <C m w , (56) 

where we have defined the natural dimensionless self coupling constant of pions by 

In this region the temperature is far smaller than the mass of massive modes (the P's and half of the Q's), but is far 
larger than the mass of the almost Goldstone modes (the other half of the Q's). We also have the intrinsic applicability 
limit of the low energy effective Lagrangian that other hadronic modes are not excited: i.e., <C Aqcd 3 - 
We expand the Laudau-Ginzburg free energy (normalized by that at /1 = T = 0) up to the quartic order in a, 

—4nF 2 M 2 =C0 + C2C? + C4 " 4 + ° (a6) - (58) 
Using the expression (|50(l , the coefficient of the quadratic term is given explicitly by 



I 2n 2 1 

C2 = -X + 



-2nM A + 4(1 + 2n)^M + 2(-nM 2 + 4(1 + n)^ 2 )Tlog(l 



2 M 2 32nF 2 M 2 

n n, , M—2u n 0l M + 2li . 

+ (-M 2 + 4fj,M + 4^ 2 )nTlog(l - e —) + (-M 2 - 4fiM + 4/x 2 )nTlog(l - e ~) 

AnfiT 2 ( Li 2 (e ~) - Li 2 (e r*)) . (59) 
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Similarly, from l|50[l the coefficient of the quartic term is given explicitly by 



c 4 = 777 



4 



10nM 3 - 56(1 + 2n)Mn 2 + 48(2 + in) — 



24 ZM 2 768ttF 2 M 2 
3n(Af 2 - AM a - An 2 ) 2 3n(M 2 + 4M/x - 4/x 2 ) 2 6 (nM 4 - 8(1 + n)M 2 n 2 + 32(1 + n)//) 
(e^ 1 - l)M (e^^ - l)M {e~ - 1)M 

+2n(Af 2 - 4A/> - 40^ 2 )Tlog(l - e"^^) + 2n(M 2 + AMu - 40/i 2 )Tlog(l - e"^^) (60) 

+4 (nM 2 ~ 16 l + nW 2 )Tlog (l-e - ^ ) + 8nnT 2 ( Li 2 (e - Li 2 e ^)) + 12rut 3 / — coth - 

V / / M-2n X 2 

The separation between the tree level part and the one-loop part is apparent in these expressions for c 2 and C4. We 
stress that no approximation has been used so far, except that (|58J) is an expansion for small values of the order 
parameter a. 

In the parameter region (|56|l we may obtain more compact approximate expressions for c 2 and C4 by the use of the 
expansion formula of polylogarithms with positive integer indices, 

luo = ^ (n) + 7 ", l0S£ (- £ )»- 1 + E c± ¥ ^(^) k - 

» ' fc>0 

and the formula 



f , x 2 , , n . 0N n /e" A 

I -coth- = - + logA + O(e0) + O^— 



(62) 



valid for < e -C 1 <C A. We define the dimensionless deviation parameter 



5 = 1-^^-2/2 (63) 

which is of the same infinitesimal order as the coupling g. Then we find the following simple approximate formulas 
for c 2 and C4: 

c 2 - -^(l-^log^), (64) 
1 bn T 

C4 °"-M +9 24MS- (65) 

In the Landau- Ginzburg approach, the critical line is characterized by the condition c 2 = 0, and the tricritical point 
is characterized by the conditions c 2 = C4 = 26] . Expressing dimensionful quantities in the unit of the renormalized 
pion mass m n , rather than the bare pion mass Af, the critical conditions given by the approximated coefficients 1651) 
read 

critical line: /1 = ngT log , (66) 

2 m v — 2/i + gm v 

tricritical point: Mtri = ^ (l + , _ \_) , ^ = - ({ . _ . (67) 

We have confirmed numerically that these expressions approximate very well the critical line and tricritical point 
derived numerically from the vanishing of the full expressions for c 2 and C4 in I|59I60|I . for the parameter region of 
concern. Accordingly our approximation is indeed self-consistent in that the critical region is within the region (|56|l . 



C. The critical line and a weakly interacting Bose gas 

The critical line l|66f) has an interpretation in terms of a weakly interacting Bose gas, as in the 3+1 dimensional case 
studied in [bH Il5| , where the diquark phase was argued to be a Bose condensate. The fact that the same structure 
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0.08 




0.4995 0.5 0.5005 0.501 

\x/m % 

FIG. 1: Plot of the phase diagram for two flavors (n = 1) and coupling g = 0.002 which corresponds to M/F 2 = 0.2. The 
solid line is the critical line in 16611 . and the x marks the tricritical point 167H . The dotted line is an approximation based on 
expanding the effective potential to sixth order in a and demanding it to have the form Q,(a) = ce,a 2 {a 2 — const.) 2 + const. 
Note that the scale has been magnified to show just the vicinity of the critical line. 



appears in our 2+1 dimensional model is quite interesting, as Bose condensation is rather different in two spatial 
dimensions. So we show here explicitly how the usual Bose gas results are consistent with the Landau-Ginzburg 
effective Lagrangian approach of this paper. 

Recall that for a free relativistic Bose gas in d spatial dimensions, the charge density is given by [2j| 



P = 



[(f) 

n (d+l)/2 



rpd 



M fi 



9d{ T ' M 9d V T ' M 



fi 



The first term corresponds to the particles and the second to the antiparticles. For d = 2 [2£ 



.92 



T M 



,(/i-M)/T 



(68) 



(69) 



In the 2 + 1 dimensional two-color QCD system discussed in this paper, the Q excitations have baryon charge 2, so 
we should replace the chemical potential in Ij68l69(l by fi — > 2^s, and the charge density in l|f)%| by p — > ub/2. At low 
temperature the lightest modes, the Q modes, dominate. There are n 2 such modes, where we recall that the number 
of flavors is Nf = 2n. Furthermore, at low temperature only the particles contribute and so we expect 



,T 2 



n B 



2n z — g 2 — 
2-k \T 



M 2^ B 



M 



(70) 



In the Landau-Ginzburg approach discussed in the previous section, the critical line separating the normal and 
diquark phases is determined by the vanishing of the coefficient C2 of the term quadratic in the order parameter a in 
the expansion of the free energy l|58() . At one-loop, the vanishing of C2 in (|59|l near the critical point fiB ~ m^/2 can 
be written as [here we use slightly less restrictive approximations than were used in obtaining (|66|) ] 



Ms 



2 



nT 
32nF 2 



Li 2 ( C (2MB-M)/T 



3 (2 MB -M)/T 



12 



Combining these two results (|70|l and (|71(l we see that 

n B = WN f F 2 (/i fl - ^) . (72) 
In the effective Lagrangian approach the renormalized vacuum energy is [3ll | 

£^ = n B — + —— + .... (73) 
where the coefficients are fixed by the symmetries. Since the baryon number chemical potential is 

we see that the relation (|72H follows directly. This demonstrates the consistency of the Landau-Ginzburg effective 
Lagrangian approach with the standard Bose gas results, and shows that the critical line separating the normal and 
diquark phases describes a free Bose gas in two spatial dimensions. 



V. CONCLUSION AND DISCUSSION 



In this paper we have used the low energy effective field theory method at finite temperature and finite baryon 
density to investigate the phase structure of three dimensional parity invariant SU(2) QCD with fundamental quarks. 
We computed the one-loop effective potential at both T = and finite T, at nonvanishing baryon number chemical 
potential. At T — 0, the tree-level critical point /i = M/2 between the normal and diquark phase receives a finite 
shift at one-loop, but this shift corresponds exactly to the mass renormalization, so that at one-loop the T = critical 
point is fi — rriir/2, where is the one-loop renormalized pion mass. At nonzero temperature we found a simple 
expression (|50f) for the effective potential, suitable for a Landau-Ginzburg analysis in terms of the order parameter a 
which describes the rotation of the tree-level vacuum away from the normal phase vacuum alignment. The subsequent 
Landau-Ginzburg analysis of this free energy revealed a critical line in the (/i, T) plane, given approximately by Ijtjtjfl . 
This line separates the normal and diquark phases, and ends in a tricritical point whose location is given approximately 
by Ht)7f) . This gene ral structure, with a critical line ending in a tricritical point, is qualitatively similar to that found 
in the 4D case [l5| , except that the T log T temperature dependence for the critical line in l|66|l becomes a power law 
- T 3 / 2 in the 4D case. 

However, the details of the 3D phase diagram are much less clear than in the 4D case, due to the presence of infrared 
divergences, which become more pronounced at higher loops, where they will appear as powers. In a restricted region 
of the relevant coupling and parameters (such as g, T/M and p), the one-loop effective potential does describe the 
correct physics. The one-loop calculation indicates that some transition happens in the vicinity of the critical line 
(|66fl . and the one-loop Landau-Ginzburg analysis suggests that this is a line of second-order phase transitions, as is 
the case in 4D [l^. However, in 3D there are infrared divergences as p, tends to zero (fi — > M/2), which become 
more significant at higher loops. This suggests that a nonperturbative analysis may be necessary to identify the true 
nature of the phase transition. There is another reason to think that in the 3D case the critical line may not be one of 
second-order phase transitions: in two spatial dimensions, only truly massless modes Bose condense - massive modes, 
no matter how light, do not Bose condense [2||. If the Q modes are truly massless even at nonzero temperature, we 
would expect to be able to dimensionally reduce the system to two dimensions, in which case the Mermin-Wagner- 
Coleman theorem tells us that a continuous symmetry cannot be broken (in the absence of long-range interactions). 
What then is the fate of the diquark phase at nonzero T? Perhaps nonperturbative effects generate an exponentially 
small, ~ e~ F l T (rather than actually zero), mass for the modes (which at tree level are exactly massless), and the 
critical line is actually a cross-over between a normal and diquark phase, until the tricritical point is reached, after 
which it presumably becomes a line of first-order transitions. We mention that such a nonperturbative exponentially 
small mass has been observed in the thermodynamics of the large N O(N) sigma model in 3D |30|. This would mean 
that the standard perturbative effective potential approach to two-color QCD in 3D is good for describing the physics 
on some portion of the critical line where it is a noninteracting Bose gas, but away from this region the perturbative 
analysis is insufficient. It would be very interesting to investigate the relevant nonperturbative effects, perhaps from 
a large Np analysis, or from a lattice study of 3D two-color QCD. 
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APPENDIX A: ONE-LOOP RENORMALIZATION 



In this Appendix we perform the one-loop renormalization of the pion mass and wave function. As the theory is 
renormalizable at each order of the chiral perturbation, renormalization constants are common for both normal and 
diquark phase. Thus we shall concentrate on the simpler case of the normal phase, and set ^4 = as mentioned in 
section IV. A. We substitute (|3"9"jl at a — into the tree-level term of l(2"9"|) and expand in P and Q fields up to the 
quartic order. We write 



Pab = ^+i^, Qab = ^+i^, 0,6 = 1,. 



Without loss of generality we choose to compute ( <pfi {p)4>\x (~p) ) ■ The terms that contribute to this Green's function 



(Al) 



4 n 



= EE ^H 2 + m2 )^ 



a—l a, 6—1 

M 2 ' 



96F 2 
1 



*)\2 



96F 2 



(1)n2 



E^V+EEfwii 

a— 2 a — l a—l 

4 n 

4£ [d^f + EE ((^ii } ) 2 + (M: } ) 2 ) 

a— 2 a—l a—l 

4 4 n 

^E^f+EE^i) 2 



There exist vertices of the form d v <j>^ <j>^ d v <p^ , a 7^ 1, but they do not contribute to 
Using these vertices, we have 

1 



{p 2 + M 2 ) 2 
1 



p z + M 



M 2 
96F 2 



96F 2 
1 



(p 2 + M 2 ) 2 
1 

(p 2 + M 2 ) 2 



2.1(12 + 8(n-l)) 

' 2 / 2n + 1 
P 1 



4-3 



d d q 



d d q 



1 



M 2 

+ 7^2-1 3 + 8 n-1 



(2?r) d g 2 + Af 2 ' 48F 2 

1 ■2-l(12 + 8(n-l))(V) 



(2ir) d q 2 + M 2 96F 2 ' 



(2?r) d g 2 + M 2 
d d q 1 



(27r) d g 2 + M 2 



d d q 



1 



1 



12F 2 J (2ir) d q 2 + M 2 
2n+lT(l-d/2 ) ATd _ 
12F 2 {4ir) d / 2 



M 1 



n — 1 



6F 2 7 (2tt)<V + M 2 



n - 1 r(l - d/2) , o 



6F 2 (47r) d / 2 



We have discarded a term proportional to the divergent integral J i ^ in the dimensional regularization 



The wave function renormalization should absorb the factor multiplying p 2 , 

<P{P) = <t>(j>)r 



x 2n + ir(l-d/2) Jlfd _ a 



24F 2 (47r) d / 2 



Accordingly, 



At p = 0, 



That is, 



</»H (P)ren^u (~p)r 



1 



(p 2 + M 2 ) 2 



p z + M z 1 - 



1 r(i-d/2) , 



tfl ) (0)ren^ff(0)ren) = 



M 2 



1 - 



4jF 2 ( 47r )d/2 

r(i-d/2)M d - 2 ' 



4 ( 47r )d/2 F 2 

* ~ ^ + 8(47r)*/ 2 F 2 J ' 



,2 ■ 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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Note that the n-dependent part of the coefficients of p 2 and M 2 in i|A3|l are identical, leading to the n-independent 
mass renormalization. 
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